Exploiting boundary states of imperfect spin chains for high-fidelity state transfer 
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We study transfer of a quantum state through XX spin chains with static imperfections. We 
combine the two standard approaches for state transfer based on (i) modulated couplings between 
neighboring spins throughout the spin chain and (ii) weak coupling of the outermost spins to an 
unmodulated spin chain. The combined approach allows us to design spin chains with modulated 
couplings and localized boundary states, permitting high-fidelity state transfer in the presence of 
random static imperfections of the couplings. The modulated couplings are explicitly obtained from 
an exact algorithm using the close relation between tridiagonal matrices and orthogonal polynomials 
[Linear Algebr. Appl. 21, 245 (1978)]. The implemented algorithm and a graphical user interface for 
constructing spin chains with boundary states (spinGUIn) are provided as Supplemental Material. 

PACS numbers: 03.67.Hk, 03.67.Lx, 75.10.Pq 
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I. INTRODUCTION 

Spin chains have attracted much attention in recent 
years (for reviews see [H, because of their ability ei- 
ther to act as quantum communication channels |3|-|5| 
or to generate highly entangled states for quantum com- 
putation The use of spin chains for both of these 
tasks has been considered in the context of various physi- 
cal systems: Implementations of spin chains may connect 
nitrogen- vacancy registers in diamond @ or entangle in- 
ternal states of an array of ultracold atoms confined to 
an optical lattice Arrays of capacitively coupled flux 
qubits have also been shown to be suited for quantum 
state transfer [loj . 

In particular, spin chains of XX type [ll| can be 
used as a quantum channel, i.e., a quantum state (qubit) 
placed at one end of the chain can be perfectly trans- 
ferred to the other end as a result of the coherent time 
evolution. Two different strategies have been suggested 
to achieve perfect state transfer (PST) in this context. 
The first approach relies on modulated couplings between 
neighboring spins throughout the spin chain 0, Q . More 
precisely, the quality of the state transfer depends on the 
energy spectrum of the spin chain and suitable couplings 
for PST are obtained by solving an inverse eigenvalue 
problem (lEP). The prototypical couplings between spins 
for PST correspond to a linear spectrum Q, also con- 
sidered relevant to the dynamics of electrons in A^-level 
systems [l^ - [T5j . 

The second approach is based on weak couplings of 
the outermost spins to the rest of the otherwise unmod- 
ulated spin chain js], ITgI [TtI . In this weak-coupling limit 
the dynamics of the spin chain reduces to an effective 
two- or three-level system consisting of states localized 
at the boundaries of the spin chain (boundary states) 
and state transfer is a consequence of Rabi-type oscil- 
lations [l^, [l^. PST is only achieved in the limit of 
vanishing couplings of the outermost spins of the chain. 
The advantage of this approach is that transfer does not 



depend on details of the spin chain, i.e., no modulation 
of inner spin couplings is required. 

The primary aim of this paper is to combine these two 
strategies. We construct spin chains with modulated cou- 
plings; however, these are chosen such that state transfer 
takes place mainly through boundary states. The result- 
ing spin chains therefore have two qualities: they permit 
PST through a perfectly engineered chain, and their dy- 
namics, involving only a few states, is robust to small 
variations of the spin couplings. To construct these spin 
chains we start from the energy spectrum and add nearly 
zero eigenvalues, which under specific conditions entails 
boundary states. 

The couplings of the spin chains are explicitly obtained 
by solving the lEP with a numerically stable algorithm 
developed by de Boor and Golub [l^ ■ We describe the al- 
gorithm adapted to the specific problem of PST through 
spin chains and provide its implementation as matlab 
code. In addition, we supply a graphical user interface 
for designing spin chains with boundary states, called 
spinGUIn, as Supplemental Material jT9| . 

As a proof of principle, we compare transfer fidelities 
of different spin chains with static imperfections of the 
spin couplings, similar to the analysis in Refs. fl7l[20l[2]| . 
The imperfections are modeled as relative fluctuations of 
the couplings drawn from a uniform distribution. The re- 
sulting distributions of the transfer fidelities reveal that 
chains with boundary states are more resilient to im- 
perfections. This is reflected in more instances of high- 
fidelity transfer through the spin chain. 

The structure of this article is as follows: In Sec. II 
we recall the basics of state transfer through spin chains. 
In Sees. Ill and IV we give necessary details of the two 
strategies for achieving state transfer using modulated 
couplings and weakly coupled end spins. Furthermore, 
we describe the algorithm by de Boor and Golub, which 
may be applied to the modulated coupling strategy. In 
Sec. V we present the combined approach and apply it to 
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specific examples, including an analysis of their transfer 
fidelity in the presence of imperfect couplings. We end 
with the conclusions in Sec. VI. 



II. STATE TRANSFER 

To start with, we recall the basics of state trans- 
fer through a onc-dimcnsional XX spin chain with the 
Hamiltonian 

N-l N 



Here, bj are the spatially dependent spin couplings be- 
tween neighboring sites, aj are local external fields, and 
a^,a'^,a^ are the Pauli matrices. It is convenient to 
map the spin Hamiltonian Hs to a onc-dimcnsional 
fermionic hopping model using the Jordan- Wigner trans- 
formation which yields the equivalent Hamiltonian 
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The operators Cj {cj) create (annihilate) a fermion at 
site j and obey the usual anticommutation relations. 
Since the Hamiltonian Hp commutes with the total num- 
ber operator n = X^jLi^jcj, the Hilbert space can be 
decomposed into subspaces Hn corresponding to differ- 
ent total fermion numbers n. For transferring a sin- 
gle qubit we restrict our considerations to the subspace 
y.o(B'Hi- The subspaces Ho and Hi are spanned by the 
vacuum state |vac) and the single-fermion Fock states 
\j) = cjjvac), respectively. 

For state transfer the spin chain is initialized in the 
vacuum state. The qubit is written into the first spin 
of the chain at time t = 0, so that the state of the 
spin chain is IV'(O)) = co|vac) + ci|l). Subsequently, 
the qubit is transferred under the coherent evolution 
U{t) = e~'^^* to the spin at site N after the transfer 
time t = T, where it can be read out. In the ideal case 
we have \^{t)) = ?7(r)|?/;(0)) = co|vacHci|A^), i.e., PST 
is achieved. The unitary operator U{t) and therefore 
the Hamiltonian Hp have to fulfill certain conditions to 
ensure PST. Since the vacuum state has a trivial time 
evolution, the conditions only apply to the Hamiltonian 
for the subspace Hi, which in the single-fermion basis 
{\j)} is given by the tridiagonal matrix 
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III. MODULATED COUPLINGS 

The approach based on modulated couplings between 
neighboring spins relies on two conditions of if at (for de- 
tails see First, the matrix Hn has to be symmet- 
ric along the antidiagonal, i.e., the entries of H^ fulfill 



the condition aAr^j+i 



and b 



N- 



bi. The matrix 



Hpf, being symmetric along the antidiagonal, is said to 
be persymmetric. As a result of the reflection symme- 
try, the eigenvectors |Afc) of the matrix H^ have definite 
parities. Moreover, if the eigenvalues Xk are in increas- 
ing order then the eigenvectors \Xk) change parity alter- 
natively, i.e., the mirror-inverted eigenstates |Afe) satisfy 
the relation \Xk) — (— 1)'^|Aa;) upon assuming that even 
(odd) k label even (odd) eigenstates |Afe). 

Second, for PST the eigenvalues A^ have to fulfill the 
condition 



= (-1)' 



(1) 



for a constant transfer time r and phase $. In fact, 
if the initial single-fermion state |<i2(0)) G Hi is ex- 
panded in terms of eigenvectors as |</3(0)) X^fc'^fel'^'^) 
with constant coefficients c^, then the state at time t 
is given by |</3(t)) = X^fe CfeC~''^''^|Afc). On the other 
hand, since \(p{t)) is the mirror-inverted state of |(y9(0)), 
have |v3(r)) = e"^Y.k^k\Xk) = 
A comparison between the two ex- 



by assumption we 
e'^E^Cki-mX,). 
pressions for |</j(t)) then indeed yields the condition in 
Eq. ([T]). The generalization by global phase factor c'* 
can be made since the phase $ can be compensated for, 
e.g., by applying a constant external field aj ^ — $/r for 
all J. 

Thus, finding the Hamiltonian Hn for PST reduces to 
an lEP, namely, calculating the couplings aj and bj for a 
given sequence of eigenvalues A^ that fulfill the condition 
in Eq. ([T]) for a fixed r. A convenient choice is to set the 
transfer time r to the fixed value tt so that the cigenen- 
ergies Afc take integer values. Other transfer times are 
obtained by rescaling the spectrum by an overall energy 
scale. 



Solving the lEP 

We now describe the algorithm developed by de Boor 
and Golub [3 for solving the lEP in the case where 
Hn is persymmetric and all eigenvalues A^ are distinct. 
Solving the lEP based on continued fractions has been 
suggested recently in Ref. [22| for achieving PST, and 
similarly in the context of electric circuit theory |23| . We 
chose the algorithm by de Boor and Golub. which was 
in part motivated by Ref. [13] , because of its clarity and 
straightforward numerical implementation. 

We start with basic definitions related to orthogonal 
polynomials and tridiagonal matrices. We denote by Hj 
the left principal submatrix, which is formed by deleting 
the last N — j rows and columns of Hn- Further, we 
introduce the polynomials Pj{x) ~ dct{x — Hj), with j = 
1, . . . , A^, and define po = 1 and p-i = 0. Clearly Pj{x) 
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are the characteristic polynomials of the matrices Hj , and 
in particular, PAr(Afe) = for the eigenvalues A^. It then 
follows directly from Laplace's formula for the expansion 
of determinants that the polynomials Pj{x) satisfy the 
three-term recursion relation 

Pj = (x - aj)pj-i - &j_iPj-2 ■ (2) 

Next we introduce the discrete scalar product ( , ) 
defined as 



(3) 



for a ny p olynomials / and g up to degree N . As shown in 
Ref. [1^ , the polynomials Pj are orthogonal with respect 
to the scalar product in Eq. ([3]), i.e., {pi,Pj) = for 
i ^ 2^ provided that the spectrum-dependent weights 
are defined by lUfc = |dp7v(a;)/dx| ^ evaluated at a; = 
Afc. Using the expression pn{x) = Y\^.{x — Xk) for the 
characteristic polynomial, one finds the explicit form for 
the weights: 



-Ap) 



(4) 



The orthogonality of the polynomials pj and the recur- 
rence relation make it possible to express the coefficients 
Qj and bj solely in terms of pj and Pj-i- By taking the 



scalar product with pj- 
obtains 



on both sides of Eq. ([2]) one 



{xPj-l,Pj^l) 



(5) 



Similarly, taking the scalar product with Pj-2 and pj 
on both sides of Eq. ([2|) yields b'j_i{pj-2,Pj-2) = 
{xPj-i,Pj-2) and {pj,Pj) = {pj,xpj_i), respectively. 
Hence using the property of the scalar product {xf, g) = 
(/, xg), one finds 



lb. 



\\P]-i 



(6) 



with thenorm||/|| = /(7J). 

The algorithm is based on the key observation that the 
polynomials Pj and coefficients aj and bj can be deter- 
mined recursively starting with the polynomials p_i = 
and po = 1- The required weights w^, which specify 
scalar product in Eq. ([3]) , are readily calculated from the 
eigenvalues Xk- Thus the algorithm for solving the lEP 
consists of the following steps: 

1. Calculate the weights Wk for the given Afc from 
Eq. ([1]). Subsequently repeat steps 2-4 for increas- 
ing j starting with j = 1. 

2. Calculate the coefficient Uj from Eq. ([5|). 

3. Find the values oipj{x) a,t x — Xk from Eq. ([2]). 

4. Calculate the coeflacient bj from Eq. 



For odd N the steps have to be repeated up to j = {N + 
l)/2 and for even N up to j = N/2. 

As noted by de Boor and Golub, the algorithm pro- 
vides a stable means of computing the entries of H^. 
However, the maximum length of the spin chain TV is lim- 
ited by floating point under- or overflow in the weights or 
the values of the polynomials. To ameliorate this prob- 
lem during execution of the algorithm it is advisable to 
scale the spectrum {A^} into the interval [—1, 1] with ap- 
propriate rescaling of the resulting coefficients aj and bj . 
This prevents the weights Wk ^ 1/AA^~^ (AA being the 
typical difference between any two of the Xk) from be- 
coming smaller than the typical floating point precision. 
Moreover, the polynomial terms Pj{Xk) in the calculation 
are potentially large without scaling since they take val- 
ues of the order of X^^^+^^Z^ for j = {N + l)/2. With 
this proviso the algorithm yields accurate results for spin 
chains with lengths up to a few hundred spins, as found 
by testing against the exact solutions for the linear and 
cosine spectrum (defined in Sec. V). 

Note that the algorithm does not involve approxima- 
tions and can be used to obtain exact analytic results 
for the coefficients aj and bj. For instance, a useful ana- 
lytic result following from Eq. ([Sj and based on symme- 
try arguments is that all aj vanish if the eigenvalues Afc 
are symmetrically distributed around zero (cf. Lemma 5 
in (H); the converse is also true [l5| . 



IV. WEAKLY COUPLED END SPINS 

The second approach relies on weak couplings of the 
outermost spins (end spins) to the rest of the chain, 
i.e., the couplings bi and b^-i are considerably smaller 
than the otherwise arbitrary spin couplings bj with j = 
2, . . . , — 2. To analyze this type of spin chain we can 
therefore use perturbation theory in the couplings 6i, b^- 
Accordingly, we partition the Hamiltonian into two parts 
Hn ^ Hq + V with 



/ 6i 
bi 



V 











V 



bN-i 




and Ho = Hiq — V , where for simplicity we assumed that 
the local external fields aj vanish. 

As customary, we introduce the eigenvalues ^fc and 
eigenstates \^k) of Hq, i.e., i/olCfe) = 6I6>- Before ap- 
plying perturbation theory we make the following obser- 
vations: First, the eigenvalues £,k of Hq are symmetri- 
cally distributed around zero. Therefore, if the number 
of spins N is odd then the Hamiltonian Hq has an eigen- 
state l^o) with = 0, i.e., a zero mode. This mode has 
the property that all odd components in the position ba- 
sis {|j)} vanish identically and that the even components 
have alternating signs. Second, Hq has two additional 
zero modes regardless of A^ which are localized at sites 
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1 and N, i.e., the two boundary states = |1) and 
I^Tv) = \N) are zero modes of Hq. In sum, the = 
subspace is spanned by the states |^o), l^i), \S,n) for odd 
N and by the states \^n) for even N. 

We determine the evolution of the state \ip{t)) initially 
in state l£i) = |1) by using time-dependent perturbation 
theory [23|- To this end \(p{t)) is expanded in the basis 
as 

k 

where Ck{t) are time-dependent coefficients with ini- 
tial conditions Cfc(O) = Sik- Inserting \(p{t)) into the 
Schrodinger equation yields 

with the matrix elements Vpk = (^pl^l^fc)- The cou- 
pled equations for the Ck{t) can be solved approximately 
by separation of time scales [l^. As a first step, the 
equations for the coefficient Ck{t) with k ^ 0,1, are 
solved under the approximation that the slowly vary- 
ing Co(t), ci(i) and c^ify are constant, which yields 
Ckit) « (1 - e'«'=*)[14ici(t) + T4.jvCjv(i)]/&- After in- 
serting these approximate solutions into the equations 
for co(t), ci(t), CAr(t) and neglecting all fast oscillating 
terms 6"'^*=* we obtain 

^7)7 ^" ~ P'oi ^ Von ] { CO ] , (7) 
\ / V 5^ An j \ cN j 

where the detunings and the frequency Q. are given 
by 

k k 

Because of the symmetry of the problem the detunings 
Ai vanish identically [l5[. We discuss the dynamics of 
the state \(p{t)) resulting from Eq. (O for the case Vbi = 
Von = vj^f^ and h\ = hN-i ~ b. 

If N is odd then the system has a zero mode |^o) 
and the dominant contributions in the limit of weak cou- 
plings 6 — > come from the matrix elements v ^ b since 
the second-order frequencies scale as ~ The state 
evolves according to 

and thus is transferred from site 1 to A^ after the time 
T = Ti/v ^ \/b. The (unnormalized) eigenstates of the 
matrix in Eq. ([7]) are given by |1) — |A^) and |1) ± |^o) + 
|A^) with eigenvalues and ±i/, respectively. Thus the 
perturbation V leads to strong mixing of the boundary 
states |A^) with the zero mode |fo) and lifts their 
degeneracy with an energy splitting proportional to b. 
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FIG. 1. (Color online) Couplings bj between neighboring 
spins for perfect state transfer before (full dots •) and af- 
ter (empty dots o) adding boundary states to the spin chain. 

(a) The outermost couplings b\, bN of the spin chain with lin- 
ear spectrum are reduced when boundary states are added. 

(b) The couplings of the spin chain with inverted quadratic 
spectrum are approximately mirrored along the horizontal 
axis, (c) Components of the eigenstates \Xk) in the position 
basis corresponding to the modified couplings in (b). 
The boundary states |A-i), |Ao) and |Ai) can be clearly seen. 
The parameters are set to A*' = 31, A = 7 and the bj are 
scaled by 1/25 and 1/100 in (a) and (b), respectively. 



If A' is even then the zero mode |^q) is absent and 
Eq. ([7]) describes on-rcsonance Rabi oscillations between 
the states |1) and |A) at the Rabi frequency Vl. Accord- 
ingly wc have 

|^(t))=cos(^^)|l)-isin(^^) |A) 

and state transfer takes place after the time r = 7r/f2 ~ 
1/6^. The eigenstates under the effect of the perturbation 
V are |1) ± |A^) with eigenvalues ztfi, and thus wc have 
an energy splitting proportional to b^. 

The perturbative approach is valid provided that 
i^, n <C ^min, where ^min IS the non-zero eigenvalue with 
the smallest magnitude. In this regime the contributions 
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from high-frequency modes with > ^nun average out 
on the time scale of the state transfer. In particular, 
PST is achieved only in the limit of vanishing couplings 
61 and b]\[-i, which is the main drawback of this ap- 
proach. On the other hand, in this limit PST is possible 
for arbitrary configurations of the inner couplings bj with 
j-2,...,iV-2. 



V. COMBINED APPROACH 

We combine the two approaches in order to endow 
modulated spin chains with boundary states. State 
transfer then takes place mainly through the boundary 
states, making it more robust to imperfections, and yet 
is perfect even for finite couplings to the end spins. As 
seen previously, nearly zero eigenvalues in the spectrum 
of the spin chain, i.e., eigenvalues significantly smaller 
in magnitude than any other eigenvalue, are a signature 
of boundary states. The crucial question is if the con- 
verse is also true, i.e., if adding nearly zero eigenvalues 
is sufficient to introduce boundary states. The answer is 
no — the presence or absence of boundary states depends 
on the entire spectrum of the chain. 

In order to show this we focus on spin chains with odd 
A'^; however, the arguments are similar for even N . The 
condition for state transfer through boundary states is 
^ ^niin, or equivalently b <C Cmin- To obtain a con- 
dition only on the spectrum we make two observations: 
First, we find from the algorithm by de Boor and Golub 
that b^ = X]fe ^fe'^fc/ Sfc '"'fcj which is the weighted vari- 
ance of the spectrum. Second, we notice that ^min is 
identical to Amin provided nearly zero eigenvalues are ex- 
cluded from the minimum. Therefore the condition for 
state transfer through boundary states in terms of the 
eigenvalues Xk and the weights Wk = | Hp^fel-^fc ~ -^p)!"^ 
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FIG. 2. (Color online) Tlie histogram and best-fit /? distri- 
bution of the overlap /i,]v for the same spin chains as in 
Fig. 1 with randomly perturbed couplings bj. (a) Adding 
boundary states to the chain with linear spectrum results 
in more instances of high-fidelity transport and a broader 
distribution of fi,N- (b) The average transport fidelity of 
the chain with inverted quadratic spectrum is significantly 
improved by adding boundary states. The best-fit /3 dis- 
tribution corresponding to the curves labeled (l)-(4) yields 
a = {96.7,19.3,1.11,15.6}, /3 = {5.88,1.63,8.65,1.48} and 
the average ^ = {0.943, 0.922, 0.114, 0.914}. 



Even though Eq. ([8|) is the desired result we gain more 
insight by discussing it on the basis of a stepwise linear 
spectrum. The specific spectrum consist of three bands 
separated by gaps F, taking the role of Amin- The m 
nearly zero eigenvalues in the central band and the M 
eigenvalues in each of the peripheral bands have inter- 
level spacing S and A, respectively. We estimate the 
values of the weights Wk to determine the regime of S, F, 
A for which Eq. ([8]) is fulfilled. Noting that the weights 
Wk depend on the distance between the eigenvalues we 
find the scalings 



1 



A^^F'"(2F) 



M 



u™" and w^'^' are the typical weights for the eigen- 
values Xk located in the central and the peripheral bands, 
respectively. Now, Eq. (jS]) is fulfilled if contributions from 
large eigenvalues in the peripheral bands to the weighted 



variance are small, i.e., w^^'^'^ <^ w^fe'^"; or equivalently 



Thus, eigenvalues in the central band lead to boundary 
states of the spin chain in the regime (5 ^ F and A ~ F. 

To illustrate the combined approach we consider two 
examples, namely, the linear spectrum and an inverted 
quadratic spectrum. Nearly zero eigenvalues are intro- 
duced by changing the original spectrum Afe to the shifted 
spectrum Afc = Afc — sgn(Afe)C, where C < Amin is a con- 
stant. This shifting procedure is applicable to all spectra 
and leaves the length of the spin chain unchanged. 

The linear spectrum is given by Xk = Ak, where the 
index k runs from — ^{N — 1) to |(A — 1), and 1 is 
an arbitrary odd integer so that r = n/A. After shifting 
with C = A — 1 the spectrum contains two eigenval- 
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ues A±i = ±1 with boundary states |A±i). Of the cor- 
respondmg couphngs only 61, bjsi-i are significantly re- 
duced, as shown in Fig.[lja). Both spin chains, with spec- 
trum Afc and Afe, support PST and their transfer times 
are related by f = At. 

The inverted quadratic spectrum is defined by A^ = 
k(N — 1 — |fc|). The corresponding couplings bj strongly 
oscillate toward the end of the chain and are approxi- 
mately constant in the center. The shifted spectrum with 
C = N — 3 has two eigenvalues A±i = ±1 with boundary 
states |A-|-i). Both spin chains, with spectrum A^, and 
Afe, support PST with identical transfer times. The cou- 
plings bj and the cigenstatcs |Afe) are shown in Figs.[IJb) 
and[ljc), respectively. 

As an aside, spin chains with cosine spectrum A^, = 
2 cos [Trk/{N + 1)] and constant couplings bj = 1 do not 
support PST but can still be endowed with boundary 
states. The shifting procedure significantly reduces the 
outermost couplings 5i, bN-i and leads to slight oscilla- 
tions of the inner couplings. This represents an alterna- 
tive to the modification of the cosine spectrum suggested 
in Ref. [IBl in order to achieve PST. 



Random static imperfections 

We now turn to the performance of spin chains in 
the presence of static imperfections in the couplings 
and show that boundary states improve the transfer fi- 
delity. For this purpose we numerically evaluate the 
transfer fidelity of spin chains with randomized couplings 
6™'^ ~ + where i? is a uniformly distributed ran- 
dom variable in the interval [— r, r]. We use the over- 
lap /i,7v(t) = |(iV|e"'-^"^|f)|, taking values in the inter- 
val [0, f], to assess the performances of the chain. The 
overlap fi^N is related to the transfer fidelity F of the 
state transfer averaged over all input states on the Bloch 
sphere by F = i -I- |/i,jv + ^fi n Q- The transfer time 
T is fixed to the value of the perfectly engineered chain 
(r = 0), in which case /i,Ar(T) = 1. 

We reconsider the spin chains with linear and inverted 
quadratic spectrum and compare the distribution of the 
overlap /i,Ar(T) with and without boundary states, sam- 
pled over 10^ transfers. In addition, we fit a /? distribu- 
tion to the numerically obtained distribution of /i,Af(''')- 
The /? distribution is defined on the interval (0, 1) by the 
probability density P{x) = [B{a, 13)]-^ x'^'^il - x)^'^, 
where a and /? are two positive shape parameters and 
B{a,/3) = J^dtf^-^l -ty-^ is the /3 function. The 
average value is given by /i = a/(a + /3) and the variance 
by a2=a/3/[(a + /3)2(a + /3 + I)]. 

Figure [5] shows the distribution of the /i,Ar(T) for the 
linear and inverted quadratic spectrum, as well as their 
shifted counterpart with boundary states. In all cases 
the randomized couplings (at the level r = 0.05) result 
in reduced transfer fidelities. Both spin chains with lin- 
ear spectrum are remarkably resilient to imperfections; 
however, the corresponding distributions of the overlap 
fi,N{T) differ noticeably, as shown in Fig. [2][a). Bound- 




Time / 

FIG. 3. (Color online) The overlap fi^M as a function of time t 
for the spin chains with the inverted quadratic spectrum of 
Fig. 1(b). The irregular time dependence of /i,jv (orange 
peaks) with high-frequency oscillations and a sharp peak at 
t = TT is smoothed out by adding boundary states, resulting 
in a regular sine dependence (blue curve) with small-scale 
oscillations (visible only as the thickness of the plotted curve). 



ary states lead to a broader distribution of /i,Ar(T) and to 
substantially more instances of high-fidelity state trans- 
fer. Spin chains with boundary states are therefore ad- 
vantageous in a scenario where varying couplings are 
caused by imperfect fabrication (or tuning) and only 
high-quality chains, say with /i,Ar(r) > 0.98, are selected. 
Note that the performance of the spin chain with a lin- 
ear spectrum can be further improved by increasing the 
parameter A. 

The fidelity of chains with an inverted quadratic spec- 
trum is clearly improved by adding boundary states. The 
average fidelities and the distributions of /i,Ar(r) differ 
significantly, as shown in Fig. [2][b) . This is readily ex- 
plained by looking at the time dependence of /i,Ar(i) in 
Fig. [21 Without boundary states the overlap /i,7v(i) car- 
ries out high-frequency oscillations and is sharply peaked 
around the transfer time r = tt. As a consequence, small 
perturbations are likely to change the position of this 
peak and to reduce the overlap /i,Ar(T) considerably. In 
contrast, /i,Ar(t) varies smoothly if the time dependence 
is mainly determined by boundary states. 

Our observations made above hold true for a wide 
range of parameters, specifically, for different levels of im- 
perfection r. This is to be expected since all spin chains 
with boundary states are described by the same effective 
two- or three-level system and therefore exhibit similar 
behavior. Moreover, random imperfections of the spin 
chain, represented by i?rnd, lead to lowest-order correc- 
tions of the states |Afe) and their energies Xk proportional 
to Ep(^fcl^rnd|Ap)/(Afe - Ap) and (Afe|i/rnd|Afc), respec- 
tively. The relevant matrix elements are small for bound- 
ary states as long as the perturbation Hrnd mainly af- 
fects the inner part of the chain, and mixing of boundary 
states with high-frequency states |Ap) is suppressed by 
1/Ap. However, if the outermost couplings are affected 
by the perturbation, strong mixing of boundary states 
may occur. 
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VI. CONCLUSIONS 

In summary, we have presented a comprehensive ap- 
proach for designing spin chains suitable for high-fidehty 
state transfer. We have combined the two strategies to 
achieve PST based on modulated couplings and weakly 
coupled end spins with boundary states. This allows us 
to exploit their respective advantages, namely, PST for 
finite couplings and resilience to imperfections. We have 
shown that a large class of spin chains can be endowed 
with boundary states by modifying their energy spec- 
trum, provided it fulfills the condition in Eq. (|5]). 

We have evaluated the performance of different spin 
chains assuming that the couplings between spins are 
affected by random static imperfections. We saw that 
adding boundary states to spin chains significantly 
changes the distribution of their transfer fidelities. De- 
pending on the specific chain, either the number of high- 
fidelity transfers or the average transfer fidelity are in- 



creased. Part of this increase is explained by the smooth 
dependence of the transfer fidelity on the evolution time 
if transfer is achieved through boundary states. 

The results for single-qubit transfer can be easily ex- 
tended to registers of qubits. Adding several nearly zero 
eigenvalues to the energy spectrum results in boundary 
states localized over a few sites at both ends of the spin 
chain. As there are many possible configurations for spin 
chains supporting PST, e.g., adapted to a specific physi- 
cal system, we provide a graphical user interface for de- 
signing spin chains (spinGUIn) as Supplemental Mate- 
rial 
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SUPPLEMENTAL MATERIAL 

We here provide a graphical user interface for 
MATLAB, called SPInGUIn, that allows the user to 
generate and modify spin chains for perfect state 
transfer. The files required to run SPInGUIn are 
iepsolve.m, spinGUIn.fig and spinGUIn.p, which 
are added as supplemental material to this article 
or, alternatively, can be downloaded from the web 
page http : / / cms .uni-konstanz . de/physik/belzig/ 
research/publications. In order to start SPInGUIn 



one has to copy the files into the Matlab current folder 
and run the function file spinGUIn.p. Furthermore, 
the implemented algorithm by dc Boor and Golub 
(iepsolve.m) is given below for quick copy and paste. 

As can be seen in Fig. [5l SPInGUIn consist of four 
main panels: The upper left panel is used to define the 
spin chain either via eigenvalues or spin couplings. The 
remaining panels show (i) the time evolution of the state 
including the transfer fidelity (ii) the tunable couplings 
and (iii) the eigenvalues and eigenstates of the spin chain. 



function [a,b,mat] = iepsolve (spec) 
% 

°/„ IEPSOLVE solves an inverse eigenvalue problem: The function constructs 
°/o the persymmetric tridiagonal matrix of size N*N for a given set of N real 
"/o non-degenerate eigenvalues. 

y. 

y. [A,B,MAT] = IEPSOLVE (SPEC) 
% 

y SPEC is the vector containing the eigenvalues (no ordering required) . 

y 

y MAT is the complete persymmetric tridiagonal matrix. 

y A is the vector containing the diagonal elements of MAT. 

y B is the vector containing subdiagonal (superdiagonal) elements of MAT. 

y 

y Supplemental material to the article by M. Bruderer et al . (2011) 

y 'Exploiting boundary states of imperfect spin chains for high-fidelity 

y state transfer' 

y 

y Determine the spectrum size N and the iteration maximum M 
N = size (spec ,2) ; 
if rem(N,2) 

M = ceil(N/2) ; 

else 

M = N/2; 

end 

y Shift and rescale the spectrum into the interval [-1,1] 
mpos = (max(spec) + min(spec) ) /2 ; 
spec = spec - mpos; 
scale = max(abs (spec) ) ; 
spec = spec/scale; 

y Initialize the necessary arrays a_j , (b_j)~2, w_k and p_{j-l} 
a = zeros(l,M) ; 
bb = zerosd.M) ; 
w = zerosd.N) ; 
pi = onesd.N) ; 

y Calculate the weights w_k from the spectrum 
tmp = spec (2 : N) ; 
for k = 1:N 

w(k) = abs(prod(l ./(spec(k) - tmp))); 

tmp(k) = spec(k) ; 

end 



°L Calculate the components a_l, (b_l)~2 and the polynomial p_l 

nom_a = sum(spec .* w); 

den = sum(w) ; 

a(l) = nom_a / den; 

p = spec - a(l) ; 

nom_b = sum(w .* p.~2); 

bb(l) = nom_b / den; 

°L Calculate the remaining components a_j , (b_j)~2 and the polynomials p_j 
for j = 2:M 
p2 = pi; 

pi = p; 

nom_a = sum(spec .* w .* pl.~2); 
den = sum(w .* pl.~2); 
a(j) = nom_a / den; 

p = (spec - a(j)) .* pi - bb(j-l) .* p2; 
nom_b = sum(w .* p.~2); 
bb(j) = nom_b / den; 

end 

7o Compile and rescale the results using symmetry properties of the matrix 
a = scale * a + mpos; 
b = scale * sqrt(bb); 
if rem(N,2) 

a= [a(l:M-l) ,fliplr(a)] ; 

b= [b(l:M-l) ,fliplr(b(l:M-l))] ; 

else 

a= [a(l:M) ,fliplr(a(l:M))] ; 
b = [b(l:M-l) ,fliplr(b)] ; 

end 

mat = diag(a) + diag(b,l) + diag(b,-l); 




FIG. 4. The authors' rendition of a spinGUIn 
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FIG. 5. Screenshot of SPInGUIn: The upper left panel is used to define the spin chain either 
via eigenvalues or spin couplings. The remaining panels show (i) the time-evolution of the 
state including the transfer fidelity (ii) the tunable couplings and (iii) the eigenvalues and 
eigenstates of the spin chain. 
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